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Using generalized Nielsen identities a formal proof is given that the fermionic quasipar- 
ticle dispersion relations in a color superconductor are gauge independent. This turns out 
to involve gluonic tadpoles which are calculated to one-loop order in a two-flavor color su- 
perconductor. Regarding the appearance of gluon tadpoles, we argue that in QCD the color 
superconducting phase is automatically color neutral. 



I. INTRODUCTION 

It is well known that sufficiently cold and dense quark matter is a color superconductor 0,0). 
At asymptotic densities, when the QCD coupling constant is small, it is possible to investigate the 
properties of a color superconductor using mostly perturbative techniques based on the fundamental 
Lagrangian of QCD, because the dominant interaction leading to the formation of Cooper pairs 
is single-gluon exchange, which is attractive in the color-antitriplet channel. The correct leading 
order result for the gap has been first obtained by Son [3J, 

^=^ M exp(~) [l + OG?)], (1) 

with c = 37r 2 /\/2- This result has been confirmed, and the constant &o nas been determined 
0, 0, @, EJ > taking into account also contributions from the quark self energy 0, |{| . 

The gap equation is usually derived from the Schwinger-Dyson equation, where some approx- 
imations are made such as the neglect of vertex corrections. In [Ir1 | it has been shown that large 



gauge dependences occur in such calculations of the value of the gap for g > g c ~ 0.8. In [11 1 
it has been shown that the corresponding gauge dependent terms appear in the gap equation at 
sub-sub-leading order, though gauge dependences arise already at sub-leading order if the gap 
parameter is not evaluated on the quasiparticle mass shell. Recently, the gap equation has been 
considered in a non-local gauge where the quark self energy and vertex corrections vanish Q to 
include higher-order corrections to the prefactor, exploiting the expected gauge independence of 
the gap as a physical quantity. 

In view of all these efforts, it seems desirable to have a formal proof of gauge independence of 
the relevant quantities. In this Note we shall consider the fermionic quasiparticle dispersion rela- 
tions in a color superconductor, which are modified by the appearance of a gap which dynamically 
breaks color symmetry. Such a proof is given in section 2 of the present paper as a straightforward 
generalization of a similar proof for finite temperature QCD [131 ] . It turns out that in order to 
formulate this proof one has to allow for the appearance of gluon tadpoles in a color supercon- 
ducting phase. We demonstrate the possible existence of such tadpoles in section 3 by a one-loop 
calculation in the example of an Nf = 2 color superconductor, discussing also briefly the issue of 
color neutrality. We argue that in QCD (in contrast to NJL models) the color superconducting 
phase is automatically color neutral. 

Our notations follow mostly those of Pisarski and Rischke, Refs. |3. Il~l. Il5| . 
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II. GAUGE INDEPENDENCE OF FERMIONIC DISPERSION RELATIONS 

We consider the inverse quark propagator in the Nambu-Gor'kov formalism [ll. 1^. . 

c-i _ ( ®+ M7o + S \ 

5 ~V <*> + + (2) 

where $ ± are the gap functions, related by = 7o[^ ,+ (Q)] t 7o, and E(Q) = C^-Q^C -1 

with the charge conjugation matrix C . Flavor and fundamental color indices are suppressed in (J2J). 
This inverse propagator is the momentum space version of the second derivative of the effective 
action, 

5 2 T 

_ (3) 
8^(x)8^(y) v=-0=^=o,^g=Ag' 

where ^ = (■ifj,ip c ) T , ^ = (i/!,$ c ), and is the expectation value of Aq, which we allow to be 
non- vanishing (see section 3). 

It should be noted that the doubling of fermionic fields in terms of ^ and ^ is just a notational 
convenience here; the effective action itself should be viewed as depending only on either (t/j, tp) or 
the set * = (ip,'ip c ) T . 

The gauge dependence identity for the effective action (generalized Nielsen identity, see ap- 
pendix) follows from considering a completely arbitrary variation of the gauge fixing function. It 
can be written as 



8T = dx 



\8ip(x) 

( (5r <5r (5r \ 

dx I . , , , 8X(m (x) + - ■ . , 8X(^ )(x) + — — j—rSX?A\(x) ) , (4) 



where the various 5X are defined in l)A3|) . 

Eq. (0J) can be cast in a more compact form using the DeWitt notation for the fermions, 

ST 



5T = r A 6X* + J dx^^SXfoix), (5) 

where i = (ip(x), ijj c {x)) T , i = ( , tp(x),ip c (x)), and the comma denotes functional derivation. Taking 
the second derivative of ©, setting xj) = = A\ = 0, Aq = Aq, and using the fact that 

6T - - =0, (6) 

^,=^=A i a =0,Ag=Ag 
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we obtain a gauge dependence identity for the inverse propagator (J5J), 

ST fj = -TrfX* + V ki 8XK + J dxJ^SXfa (x). (7) 

Up to this point, our functional relations are completely general and apply also to the case 
of inhomogeneous condensates which may be realized in the so-called LOFF phase (see for a 
recent review) . In this Note we do not attempt to cover the complications this case may add to the 
question of gauge independence, but continue by assuming translational invariance. This allows us 
to introduce 5A a0 := — SXffiJx = 0) and to write © as 

dV ■-■ 

5T t3 = -Y k3 5X^ - 8X% k - q^SA*. (8) 
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Furthermore, we can transform eq. (JHJ) into momentum space, 

STfJQ) + SA^^M^ = -T kl {Q)5X%Q) - 5XK(Q)T fk {Q), (9) 



5det(r 4j ) + <Si a0 -^ det(r iJ ) ee Stat det(T^) = - det(T i5 )[5^ + <5X fc 5 ]. (11) 



where the indices i and i from now on comprise only color, flavor, Dirac and Nambu-Gor'kov 
indices. Using the fact that 

<5detM = (detM)Ti[M~ 1 8M] (10) 
for any matrix M, we obtain from eq. jUJ) 

d_ 

The left hand side of this identity is the total variation , ^| of the determinant of the inverse 
quark propagator, with the first term corresponding to the explicit variation of the gauge fixing 
function, and the second term coming from the gauge dependence of Aq. 

Since the determinant is equal to the product of the eigenvalues, eq. (|llj) implies that the 
location of the singularities of the quark propagator is gauge independent, provided the singularities 
of SX k k do not coincide with those of the quark propagator. As in |l3| one may argue that SX 
is 1PI up to a full ghost propagator, and up to gluon tadpole insertions (see Ref. 0] for the 
explicit diagrammatic structure). The singularities of the ghost propagator are not correlated 
to the singularities of the quark propagator. Their fundamental difference is in fact enhanced 
by medium effects — e.g., there exist no hard-thermal/dense- loop (HTL/HDL) |2(j corrections to 
vertex functions involving ghosts, whereas all physical degrees of freedom acquire HTL/HDL self- 
energy corrections. 

Gauge independence of the zeros of the inverse fermion propagator then follows provided that 
also the 1PI parts of 5X have no singularities coinciding with the singularities of the propagator. 
An important caveat in fact comes from massless poles in the unphysical degrees of freedom of the 
gauge boson propagator, which are typical in covariant gauges and which can give rise to spurious 
mass shell singularities as encountered in the case of hot QCD in [2l|. But, as was pointed out in 
these apparent gauge dependences are avoided if the quasiparticle mass-shell is approached 
with a general infrared cut-off such as finite volume, and this cut-off lifted only in the end, i.e., 
after the mass-shell limit has been taken. The tadpoles do not introduce any singularities. 

The determinant which appears in (jllj) is taken with respect to color, flavor, Dirac and Nambu- 
Gor'kov indices. The determinant in Nambu-Gor'kov space may be evaluated explicitly by noting 
that 

det(XU - VBV^C) (12) 

for arbitrary n x n matrices A, B, C, T> (with D invertible). Hence we obtain from (j2j) 

det(Ty) = det[($-^ 7 o + S)($+^7o + 5]) 

-($-^70 + £)$-($- ^o + S)- 1 ^]. (13) 

The inverse of the matrix of which the determinant is taken here appears, of course, in the ordinary 
quark propagator, which is obtained by inverting (|2"|). 

G+(Q) = W- ^70 + £)(#+ ^70 + S) 

-(#- M7o + ^70 + Sr^T 1 ^" ^70 + E). (14) 
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At leading order, when the quark self energy can be neglected, eqs. (fTTj) and (fT3)) imply that 
the gap function is gauge independent on the quasiparticle mass shell. It should be noted, however, 
that at higher orders only dispersion relations obtained from (|13|) . which also include the quark 
self energy S, can be expected to be gauge independent. 



III. GLUON TADPOLE AND COLOR NEUTRALITY 



As we have seen, the gauge dependence identities involve gluon tadpoles in a non-trivial 
manner. 

In a color superconductor global color symmetry is dynamically broken by the diquark conden- 
sate. This means that there is no symmetry which forbids the existence of gluon tadpoles. Let us 
consider the one-loop tadpole diagram, 



T ° = "I / ^ T ^f^KS(Q)}- (15) 



Here S(Q) is the quark propagator in the Nambu-Gor'kov basis, whose inverse is given by (sup- 
pressing color and flavor indices) 

S-^t^J^J' ) (16) 
with <I>~ = 7o( < £ + )^7o- The quark-gluon vertex Tq is given by 

f 8=( r o 8 f° 8 )- < 17 > 

with Tq = 7oT a and fg = — 7o(T a ) T . The trace in (|15|) has to be taken with respect to Dirac, 
color, flavor and Nambu-Gor'kov indices. 

The gap <3? + is a matrix in Dirac, color and flavor space. We consider as an example an Nf = 2 
color superconductor and take the following ansatz for the gap 0, |2j|: 

*7o,ij(Q) = Wi*(<l>+{V+ + - V+) + 4>-{V-_ - Vf + )), (18) 

where /, g are flavor indices and i, j are fundamental color indices. The Vs are the projection 
operators introduced in We have assumed for simplicity that the right-handed and left-handed 
gap functions are equal up to a sign 0|. 

First we evaluate the trace over Nambu-Gor'kov space which gives 

ja = -\\ ^ TT D,cj[r a g + (Q) + ng-(Q)], (19) 

where is obtained from inverting (|16|) jl~i]. 

Sf g ,ij = S fMl - + 6 i3 S j3 G±], (20) 

with Q 

^° ^-lViff A ^ (21) 

G » (<3) = £f^F A±e(q) * (22) 
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where A 1 * 1 are energy projectors, A ± (q) = ^(1 ±7o7*(f ). Evaluating the trace over flavor and color 
space we get 

r a = g(T a ) 33 J ^ I Tr D [ 7 o(G+ - CT - G+ + G )]. (23) 



Assuming that <f> ~ jj| and that 4> + has negligible four-momentum dependence in the vicinity 
of the quasiparticle pole we obtain 



T a ~ -Ag{T a ) 33 J ■ 



d^Q ( /j, — q fx — q 



*(2vr)4 Wo - (M- Q) 2 ~ I0 + I 2 Q 2 o Q) 2 . 

* ~ 2 (T a U rdqq 2 ^-q)[ 1 (24) 

TT 2 JO \ y/(fi-q) 2 + \(j) + \ 2 \H-q\J 

In order to obtain an order of magnitude estimate, we make the approximation (j) + (q) — (j)Q6{2p,—q) 
with <j)Q = const. ^J. Then the g-integration can be readily performed, with the result 

T a „ -^(T-) 33 a (cP+f In + 0((g^+f), (T a ) 33 = -5 a8 /V3. (25) 

This result is in fact of order ucj) 2 , because ln(^>/(2/i)) is of order 1/g. At this order it cannot be 
excluded that there are cancellations from higher-loop contributions which have been neglected so 
far ^4]]. In particular it is still possible that T a = 0, although this seems to be unlikely as there is 
no symmetry which forces T a to be exactly zero. 



To address the question of color neutrality 24, 25, 26, 2'3] we consider the partition function 



exp(-fi/T) = i vip exp(-S[<p\), (26) 

where (p denotes the set of all fields, and S\y>] is the QCD action (including gauge fixing terms and 
ghosts). Following the argument given in |28| it is easy to see that the system described by this 
partition function is color neutral, at least if one chooses a gauge fixing which does not involve Aq, 
for instance Coulomb gauge: 

The fields Aq appear in the action as Lagrange multipliers for the Gauss law constraint [2{|. 
Therefore in the path integral the integration over the zero- momentum modes Aq^ =0 produces 
delta functions, 5(N a ), where are the color charges. This means that only color neutral field 
configurations contribute to the partition function, q.e.d. 

The gluon tadpole obtained in (|25[) may in fact be interpreted as an effective chemical potential 
for the color number eight of the order 

gm~^n4> 2 ~ 4> 2 /(gv), 

where oc g/i is the corresponding leading-order Debye mass It may be noted that the 

chemical potential /is which has been found by requiring color neutrality in an NJL model is also 
proportional to 4> 2 [26j. We emphasize however that whereas in NJL models color neutrality has 
to be imposed as an additional condition, color neutrality is guaranteed automatically in QCD by 
the integration over the Aq zero-momentum modes. 



IV. CONCLUSIONS 



We have presented a formal proof of gauge independence for the fermionic quasiparticle disper- 
sion relations in color superconducting QCD. As long as the quark self energy can be neglected, 
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this implies gauge independence of the gap function on the quasiparticle mass shell. In general, 
however, only the singularities of the propagator, which involve also the gauge-dependent quark 
self energy, can be expected to be gauge independent. At higher orders, gauge parameter variations 
may also affect the temporal component of the gluon tadpole, which, as we have shown in a simple 
example, can be nonvanishing in a color superconductor. While a nonvanishing gluon tadpole may 
be viewed as an effective chemical potential for color charge, we have argued that color neutrality 
is an automatic consequence of QCD, if no explicit chemical potentials for color are introduced. 
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APPENDIX A: GAUGE DEPENDENCE IDENTITY FOR THE EFFECTIVE ACTION 



An arbitrary gauge theory is defined by an action Si nv [(pj which is invariant under (infinitesimal) 
gauge transformations 1 , 

5^ = DUm a - (Al) 



In order to quantize the theory it is necessary to fix the gauge freedom, for instance with a quadratic 
gauge-breaking term, 

S[<p] = S inv [y] + ±F a [<p}F a [<p]. (A2) 

By performing a certain (non-local) gauge transformation on the path-integration variable <p [l^l . 
|30j, it can be shown that an infinitesimal variation of the gauge fixing function, 5F@, induces a 
change in the effective action as follows, 

5T[ip] = -r < [0\ {D^ [<p] g a p [<p] 8F P [<p] ) [<p] =: TMSX'm, (A3) 

where Q a p\p\ is the gho st propagator in a background field <p. This gauge dependence identity is a 
generalized version of the Nielsen identity Ilia, |3jJ . 
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